Multi-component magneto-optical conductivity of multilayer graphene on SiC 
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Far-infrared diagonal and Hall conductivities of multilayer epitaxial graphene on the C-face of 
SiC were measured using magneto-optical absorption and Faraday rotation in magnetic fields up 
to 7 T and temperatures between 5 and 300 K. Multiple components are identified in the spectra, 
which include: (i) a quasi-classical cyclotron resonance (CR), originating from the highly doped 
graphene layer closest to SiC, (ii) transitions between low-index Landau levels (LLs), which stem 
from weakly doped layers and (iii) a broad optical absorption background. Electron and hole type 
LL transitions are optically distinguished and shown to coexist. An electron-hole asymmetry of the 
Fermi velocity of about 2% was found within one graphene layer, while the Fermi velocity varies 
by about 10% across the layers. The optical intensity of the LL transitions is several times smaller 
than what is theoretically expected for isolated graphene monolayers without electron-electron and 
electron-phonon interactions. 



I. INTRODUCTION 

Graphene attracts much interest due to the intrigu- 
ing physical properties and a potential for novel appli- 
cations. Epitaxially grown graphene on SiCP is partic- 
ularly promising for large-scale production. The growth 
conditions of epitaxial graphene were constantly improv- 
ing in the last years, which resulted in a macroscopic 
continuity of graphene layers and an enhanced mobility 
of charge carriers^^. Electronic properties of epitaxial 
graphene grown on Si- and C-faces of SiC, are markedly 
different!^. Interestingly, multilayer epitaxial graphene 
produced on the carbon face shows a number of electronic 
features typical of isolated monolayer graphene as re- 
vealed by infrar ed sp ectroscopy^, scanning tunneling mi- 
croscopy (STMJpSM angle-resolved photoemission spec- 
troscopy (ARPESJp) and quantum Hall effeclP. In par- 
ticular, the Landau levels in this material demonstrate 
a square-root dependence on the perpendicular magnetic 
field B and the level index n(— 0, ±1...): 



E n = E D + sga(n)^2efiv F \nB\, (1) 

where vp is the Fermi velocity and Ejj is the Dirac-point 
energy. It is generally believed that such an effective 
electronic interlayer decoupling is a result of twisted, 
non-Bernal stacking of the C-face grown graphene lay- 
ers. Significant progress in theoretical understanding of 



the influence of stackingon the electronic structure and 
Landau levels was made^3H22l. However, no complete the- 
ory able to quantitatively predict the effect of twisting on 
the band structure at an arbitrary stacking angle exists 
at the moment even for bilaycr graphene. In real sam- 
ples, many layers are present with a random rotation 
between each pair of neighbors. Moreover, the substrate 
induces a strong variation of the Dirac-point energy Ed 
with respect to the chemical potential and therefore a 
very different density and mobility of carriers in different 
layers. More experiments are needed to understand the 
complex electronic structure of this system and establish 
favorable conditions for applications. 

Infrared spectroscopy, which is a direct probe of charge 
dynamics, is well suited to study carriers in graphene. 
Optical spectra contain contributions from all graphene 
layers, including those that are not seen by surface 
probes, such as ARPES and STM. In a recent worlP^, 
we measured the rotation of polarization of light passing 
through epitaxial graphene in a magnetic field, known 
as the Faraday effect. Being an optical analogue to the 
d.c. Hall effect, the Faraday rotation reveals the sign 
of the charge carriers involved in various LL transitions. 
Moreover, optical spectroscopy allows separating carriers 
by their high-frequency dynamics, in contrast to the d.c. 
transport measurements. Here we present an extensive 
magneto-optical study of multilayer epitaxial graphene 
grown on the carbon side of SiC, by combining the Fara- 
day rotation measurements with the transmission spectra 
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and extracting both the diagonal and Hall optical con- 
ductivity. Different magneto-optical contributions are 
disentangled and studied quantitatively using a multi- 
component CR model. This approach reveals a num- 
ber of interesting properties of multiple charge carriers 
in multilayer graphene, for example, the simultaneous 
presence of both electron- and hole-like components as 
well as the existence of carriers with different Fermi ve- 
locities. It also allows us to analyze optical intensities 
of various transitions and make comparison to theoret- 
ical models. Finally, the dependence of optical spectra 
on temperature and environmental (surface) doping are 
studied. 



II. TECHNIQUES 

The rotationally stacked multilayer graphene on the 
carbon terminated side of 6H-SiC was produced by hy- 
drogen etching at 1450°C and subsequent graphitization 
in an argon atmosphere at 1650°C. The substrate has a 
thickness of 370 /mi and a surface area of 10x10 mm 2 . 
The back side of the substrate was cleaned from unde- 
sirably grown graphene using scotch tape. Using X-ray 
photoemission spectroscopy we estimated the number of 
graphene layers to be 5 ± 1 and checked that the back- 
side is graphene free. Additionally, we performed mid- 
infrared transmission microscopy over an area compara- 
ble to the total area used for the magneto-optical exper- 
iment using a spot size of 10x10 /mi 2 . We assume that 
each graphene layer absorbs 1.5% of light at 1 eV, which 
is the value expected for monolayer graphene on SiC, 
provided that the chemical potential does not exceed 0.5 
eV. This absorption value is lower than the absorption 
of 2.3% of free-standing graphene, because the refractive 
index of SiC is larger than 1. The infrared absorption 
shows a distribution peaked at 7 layers with a FWHM of 
2.7 layers. The slightly larger number of layers revealed 
by infrared microscopy can be related to an underesti- 
mation of the absorption of one layer. Notably, both 
measurements show that at least 4 layers are present at 
any point of the sample. 

We measured the Faraday rotation 9{w) and magneto- 
optical transmission T(u) in the range of photon energies 
tku between 8 and 80 meV as described in Ref. [231 The 
optical transmission was measured with respect to a bare 
SiC substrate without graphene, which underwent simi- 
lar hydrogen etching at 1450°C as the graphitized sam- 
ple. The optical spot had a diameter of 5 mm. Two series 
of measurements were performed: first by varying tem- 
perature between 5 and 300 K at a fixed magnetic field of 
3 T and then by varying magnetic field from to 7 T at a 
constant temperature of 5 K. Within each series, the op- 
tical absorption and the Faraday rotation were measured 
one after another with a shortest possible delay. During 
the experiments the sample was in a He gas flow; between 
the production and the experiments, as well as between 
the first and second series, it was stored in desiccated 



air. As will be shown later, even these precautions did 
not allow us to completely avoid surface contamination, 
which has an effect on charge dynamics in the top layers. 
Therefore, the data presented in this paper differ some- 
what from data presented in RefPS, although the same 
sample was used. 

The real parts of the diagonal, cr xx (ui), and Hall, 
a xy (ui), optical conductivities can be directly obtained 
from the magneto-optical absorption and Faraday rota- 
tion by using the general thin-film approximation, taking 
internal reflections in the substrate into account: 

l-T(w)w22 , o/.(w)Re[o- aaB (w)] (2) 
0(w) ps Z f s (u)Re[a xy (u)} (3) 

where Z ps 377 is the impedance of vacuum and f s (ui) 
is a spectrally smooth dimensionless function specific to 
the substrate^!. These relations are obtained by linear 
expansion of the exact Fresnel formulas, which is accu- 
rate for our sample with only a few graphene layers. We 
reduced the spectral resolution to 1 meV in order to sup- 
press the Fabry-Perot interference in the substrate. In 
this case: 
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where q = [(n s — l)/(n s + 1)] exp [— (2ui/c)k s d], d, 
n s (uj) and k s (co) are the thickness, the refractive index 
and extinction coefficient of the substrate. The latter 
quantities, which are independent of magnetic field, were 
determined from measurements of the absolute transmis- 
sion and reflection spectra of the bare SiC at every used 
temperature. 

The complex magneto-optical conductivity tensor: 



(oj) cr X y(u!) 
(uo) <t xx (lu) 



(5) 



is diagonal in the circular basis x ± iy, with eigenvalues: 



a± (u) = a xx (oj) ± ia xy (u) . 



(6) 



The absorption of right- and left-handed circular polar- 
izations (which are shown in the insets of FigJTJ: and [TJi) 
is described by the real parts of a + (uj) and g-JZj) respec- 
tively. This basis provides an intuitive way to present 
the magneto-optical conductivity showing individual LL 
transitions, which are active for strictly one circular po- 
larization or the otheiPffl^l In particular, for right cir- 
cularly polarized light, transitions between electron like 
LLs with n > 0, such as LLg^!, LL 1 _ > 2 etc, are active, 
provided that the chemical potential is in between the 
corresponding LLs (the symbol LL^j is used to desig- 
nate the LL transition between levels i and j). Similarly, 
for left circularly polarized light, transitions between hole 
like LLs, (LL_i^o> LL-2->-i, etc.) are excited^. 

From Eq. Q one can see that the imaginary part of 
a xy (uj) is needed to obtain the real part of the magneto- 
optical conductivity in the circular basis. Experimen- 
tally, lm[a xy (uj)] is directly related to the ellipticity of 
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the transmitted light. However, in the present experi- 
ment the error bars on the ellipticity are larger than the 
ones on the Faraday angle. Therefore, instead of directly 
applying Eq. ^ to the experimental data, we plot a± {us) 
obtained from the results of the mult i- component fitting 
discussed in Section |III A| For simplicity, in the rest of 
the paper, the symbols cr xx , cr xy , a + and cr_ will be used 
to denote the real parts of the corresponding complex 
functions, unless stated otherwise. 



III. RESULTS 
A. Magneto-optical conductivity at 3 T and 5 K 

Figure [T] shows representative conductivity spectra at 
B = 3 T and T = 5 K. The diagonal and Hall con- 
ductivities are plotted as solid lines in Fig. [T^, and [Tp. 
They arc expressed in dimcnsionless units of uq = e 2 /4ft, 
which is equal to the universal optical conductivity of 
monolayer graphene^. The rich structure of the spec- 
tra indicates the presence of multiple optical transitions. 
The absorption at low energy, shown by the arrow, cor- 
responds to quasi-classical CR coming from the highly 
doped graphene layer closest to the SiC substrate^. At 
about 60 to 70 meV, a strong peak in a xx (u>) is observed 
that matches the energy of the LLo_>i or LL_i^o tran- 
sitions at this field. The Hall conductivity displays a 
'zig-zag' shape in this spectral range, suggesting that 
there are multiple components contributing to the opti- 
cal response. In addition, a structure is present at about 
27 meV, most clearly seen in a xy (u>), resulting from the 
LLi_y2 transition. 

In order to disentangle different contributions to the 
magneto-optical spectra we used a multi-component 
model where the total conductivity is given by a sum 
of separate cyclotron resonances. In the circular basis 
the total complex magneto-optical response is given by: 



2W j 



to =F u c j + i-yj 



(7) 



and the complex diagonal and Hall optical conductivities 
are accordingly expressed by: 
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2Wj ~ w ej 



n - ( w + *7j) 2 ' 



(8) 
(9) 



where Lo c .j is the cyclotron frequency, Wj is the spec- 
tral weight and jj is the broadening of the jth compo- 
nent. Physically speaking, each component can describe 
either a quasi-classical CR or a transition between indi- 
vidual LLs separated by energy |ftaj CJ |. Notably, a xx (uj) 
does not depend on the sign of the charge carriers, unlike 
<j X y(cu), where the sign can be derived directly from the 
spectral shaped In a xx (u) each component results in 



a peak centered at |w c j|. In cr xy (uj) however, it shows 
an antisymmetric structure, where |w CJ | corresponds to 
the inflection point in the curve. The slope at this point 
coincides with the sign of the cyclotron frequency and 
reveals the polarity of the charge carriers involved in the 
transition^. In particular, a positive (negative) slope 
signals electron (hole) like carriers. 

The spectra a xx {uj) and a xy (u>) were simultaneously 
fitted using Eqs. ([8| and while allowing the parame- 
ters u> c ,j, Wj and jj to change freely. The fitting curves 
are shown in Fig. [IJi and[]J>, as black dotted lines. We 
found that a minimal model describing satisfactorily the 
spectral structures contains six components, which are 
shown separately in the same panels. 

Interestingly, it was necessary to introduce at least 
three components with transition energies between 60 
meV and 70 meV to describe the structure in this range. 
Two of the resonances are electron like, which we des- 
ignate as LL _j. l a and LLo_>i.h, and one is hole like, 
referred to as LL_i^o- Note that all these transitions 
are actually at different energies, which is clearly seen 
already from the presence of three inflection points in 
the optical Hall conductivity. A fourth component is 
at about 27 meV and corresponds to the electron like 
LLi_j.2 transition. The fifth component with a small 
value of hu c = 9 meV originates from a quasi-classical 
electron like CR. Finally, there is a component with zero 
cyclotron frequency and large scattering, which forms a 
broad absorption background present in <j xx (uj), but ab- 
sent in a xy {uj). A possible origin of this background will 
be discussed in Section ITVl 

The different contributions from electrons and holes to 
the magneto-optical conductivity are easily seen in the 
basis of circularly polarized light as shown in Fig. [TJ; 
and [TJi. The elect ron like transitions LLo_>i i(J , LLo->-i,b> 
LLi_j.2 as well as the CR show peaks in er + (w), while 
the hole like LL-i^o transition manifests itself in cr_(cj). 
Because of the small value of tu c and a relatively large 
scattering, the CR component has also a tail in ct-(uj). 
The absorption background contributes equally to cr+ (cj) 
and er_(aj). 



B. Dependence on magnetic field 

The magnetic field dependence of a xx (uj) and a xy (uj), 
measured at 5 K, is shown in Fig. [2^ and[2}}. The multi- 
component character of the spectra is present at different 
fields, although the resonances show a strong field depen- 
dence. The model described in Section UlI Al was used to 
fit cr xx (uj) and a xy (uj) simultaneously for every field. The 
fits, shown as black dotted lines, are used to calculate the 
magneto-optical conductivities in the basis of right and 
left handed circularly polarized light according to Eq. 
which are plotted in Fig. [2J; and[2ji. 

The field dependence of all fitting parameters is shown 
in Fig. |3| The energy of the CR is linear in field, as seen 
in Fig. [3k. This quasi-classical behavior is expected for 
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FIG. 1. The diagonal (a) and Hall (b) optical conductivity at 3 T and 5 K, normalized to the universal conductivity cro, and 
results of the multi-component modeling described in the text. The components in the graph are the individual contributions 
needed to create the total fit shown by the black dotted line on top of the data. The components in panel (b) have an offset 
indicated by the blue dashed line at the right side of the graph. Panels (c) and (d): the optical conductivities in the circular 
basis o + (w) and <J-(u) corresponding to the absorption of right and left circular polarized light respectively. In the insets the 
low-index LL transitions allowed in the corresponding polarizations are shown. 



highly doped graphene, where large-index LLs are close 
to the Fermi energyBSBSl j n this case the quasi-classical 
cyclotron frequency is inversely proportional to the Fermi 
energy e F : 



\e\Bvl 



(10) 



Using this relation and taking as an estimate Vp = 1.0 x 
10 6 m/s we obtain e F = 0.24 eV, which corresponds to a 



carrier concentration n 



4.2 x 10 1 



The spectral weight W of the cyclotron peak, shown in 
Fig. [3ji, is field independent within the experimental ac- 
curacy. In the absence of interactions, the Drude weight 
in a single graphene layer is related to e F by the formula: 
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This gives e F = 0.16 eV and accordingly n = 1.9 x 
10 12 cm -2 , which are significantly smaller than the values 



based on the cyclotron frequency. It is possible that in- 
teractions renormalize the Drude weight and spread the 
missing weight over a large spectral range. In our fits 
this missing weight might be 'absorbed' by the broad 
background component mentioned above. On the other 
hand, we obtain a reduced Drude weight based on the 
assumption that for the bottom layer vp = 1.0 x 10 6 
m/s. Accepting a smaller Fermi velocity (by about 
20%) would make the cyclotron resonance and the Drude 



weight match according to Eqs. (10) and (11 
The scattering rate ftry (Fig. IjjT 



otron peak 



of the cycl 

is field independent and equal to about 14 meV. Using 
the semi-classical relation: 



M = 



7 |S| 



(12) 



we find a mobility of about 2000 cm 2 /Vs, based on the ex- 
perimental values of the scattering rate and the cyclotron 
frequency. Magneto-transport measurements show that 
the carrier density and mobility in multilayer graphene 
are correlated^. The values of n and /i that we obtain 
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FIG. 2. Magneto-optical conductivity of multilayer graphene at 5 K, normalized to the universal conductivity uo, for several 
magnetic fields up to 7 T. The curves in all panels are offset as indicated by the dashed lines. Larger offsets correspond to 
higher magnetic fields. Panels (a) and (b) show the measured spectra of a xx {uS) and er xy (u)) (solid lines) and multi-component 
fits (black dotted lines). In panels (c) and (d) the model-derived o+{w) and <t_(cj) are shown. 



Transition 


VF 


[id 


AO 


LLo->-i,a 


1.02 


X 


10 B 


LLo->i,6 


1.11 


X 


10 6 


LL-i^o 


1.09 


X 


10 6 


LLi^2 


1.01 


X 


10 6 



TABLE I. Fermi velocities found for different LL transitions. 



here for the CR fall on the generic dependence found in 
Ref. [301 CR was also observed in monola yer ep itaxial 
graphene on the silicon side of silicon carbide^^H where 
similar mobility of carriers was founcP^. 

Next we discuss the field dependence of the LL tran- 
sitions. From (J+{uj) and <J-(oj) it is particularly evi- 
dent that the two electron like and one hole like reso- 
nances described above have a similar magnetic field de- 
pendence. The energies of all low-index LL transitions, 
namely LL X _>2 (Fig. |3ja,), LL ^.i, o , LL ^i !b and LL_i^ 
(Fig. [3Jd), clearly follow the square-root dependence on 
magnetic field typical of massless Dirac fermions (Eq. 
([lj), in accordance with previous observations in simi- 
lar e pitaxi al graphene samples^ and exfoliated monolayer 
flakes22E21. Fitting the field dependence of the LL tran- 
sition energies using Eq. provides the Fermi veloci- 
ties listed in Table [IJ they show a spread of about 10%. 
The difference between the Fermi velocities found from 
the electron like LLo-^i^ and hole like LL_ 1 _ ! .o transi- 
tions is only about 2%. However, the presence of three 
distinct inflection points in a xy (uj) corresponding to the 
energies of the three LL transitions, is a clear sign that 
the transition energies and therefore the Fermi velocities 
are different. 



The spectral weights of the LLi_>.2, LLo->i. a , LLo-^f, 
and LL_!__j.o transitions increase with magnetic field (Fig. 
[3ji and [3^) . This is in a qualitative agreement with an 
increase of t he nu mber of states within each LL with 
magnetic fielcf^iJ. However, a closer look at the abso- 
lute values of the spectral weights reveals a strong reduc- 
tion with respect to the theoretical expectation for ideal 



graphene monolayers, as discussed in Section IV 

The scattering rate of the LL transitions is about 3 to 
4 meV and magnetic field independent (Fig. [3^ and |3jr) . 
Although in the case of the cyclotron peak we extracted 
the mobility using the scattering rate and the cyclotron 
frequency, in the quantum limit where the transition en- 
ergy is proportional to \[B, the quasi-classical Eq. ( 12 ) 



will formally result in diverging mobility values at low 
fields. 



C. Temperature dependence 

The temperature dependence of a xx (uj) and a xy {ui) at 
a fixed field of 3 T is shown in Fig. [4^ and |4]d, respec- 
tively. In both sets of curves, the CR structure changes 
weakly with temperature. The situation is different for 
the LL transitions. In a xx (oj) the LLs are preserved at 
all temperatures, although at elevated temperatures the 
peaks are about 3 meV broader than at low tempera- 
tures as obtained from the fits shown in the same graphs. 
Their spectral weight is almost temperature independent. 
However, in a xy (uj), the spectral features corresponding 
to the LL transitions are strongly diminished at high tem- 
peratures and almost disappear at room temperature. 
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FIG. 3. Parameters of various components as a function 
of magnetic field obtained from the multi-component fits to 
a xx (ui) and a xy (uj) of Fig[2] (a) The transition energies of 
the CR (triangles) and the L Li^ 2 transition (diamonds) . The 
solid lines are fits using Eq. and Eq. |T| respectively, (b) 
The transition energies of LLo->i,o (squares), LLo_n,b (trian- 
gles) and LL_i_>o (circles) transitions. The solid lines are fits 
using Eq. |TJ. (c) Spectral weight of the LLi_,2 transition, 
(d) Spectral weight of the CR peak, (e) Spectral weights of 
the LLo-»i, a , LLo_>i,b and the LL_i_+o transitions, (f) The 
broadening of the LLi_>2 transition, (g) The broadening of 
the CR transition, (h) The broadening of LLo-»i, a , LLo-^i, 
and the —I — > transitions. Errorbars in all figures are the 
standard deviation of several fitting results. The error bars 
in panels (a) and (b) are within the symbol size. 
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FIG. 4. Magneto-optical conductivity of multilayer graphene 
at 3 T, normalized to the universal conductivity <7q, for several 
temperatures between 5 K and 300 K. The curves in all panels 
are offset as indicated by the dashed lines. Larger offsets 
correspond to higher temperatures. Panels (a) and (b) show 
the measured spectra of o xx (us) and a xy (cj) (solid lines) and 
multi-component fits (black dotted lines, only at 5, 100 and 
200 K). In panels (c) and (d) the model-derived u+(lj) and 
a-(uj) are shown. 



D. Effects of environmental doping 



The extinction of the spectral structures correspond- 
ing to the LL transitions in a xy (uj) with warming up is 
due to the simultaneous presence of electron and hole like 
LL transitions as seen in the graphs of <j+{uj) and a-(u>), 
plotted in Fig. [3J: and |4]i. At high temperatures, the 
electron- and hole- like components that have slightly dif- 
ferent Fermi velocities overlap more, due to the increased 
broadening. In <7 xy (ui) the contributions are subtractive 
and more overlap results in weaker spectral structures. 
On the other hand, in <j xx {uj), which is unsensitive to the 
sign of the charge carriers, the contributions are additive 
and the transition peak appears less affected. 



A careful inspection of the two series of measurements 
(Figs. [2] and [4]) shows that the spectra from the different 
series taken at the same experimental conditions (5 K 
and 3 T) are not precisely the same, especially close to 
the LLo-^a, LLo-^i^ and LL_i^o transition energies. 
As was mentioned above, between these measurements 
the sample was held in dry air, therefore it is reasonable 
to assume that the difference is due to the effect of the 
environmental molecular contaminatiorP^. In the case 
of multilayer graphene, the electronic properties of the 
outmost layer are most strongly modified, although one 
cannot exclude a certain effect on inner layers as well. 

Fig. [5k shows the contributions of the LLo-i-i^ and 
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FIG. 5. The effect of environmental doping on the LL transi- 
tions at B = 3 T and T = 5 K. Panel (a): the contributions of 
the LLo^i.a and LLo_>i,6 transitions to <t+(uj). Panel (b): the 
contribution of the LL_i_>o transition to <r_ (ui). In both pan- 
els, the red dash-dotted lines correspond to the measurement 
series as a function of temperature at constant magnetic field 
(Fig. 0, while the blue solid lines correspond to the mea- 
surement series as a function of magnetic field at constant 
temperature (Fig. |2J. The curves are derived from the multi- 
component fits of experimental data. 



LLo-^b transitions to ct+(cj), for both measurements at 
5 K and 3 T. The low energy LLo_yi, Q transition with Dp 
= 1.02 x 10 6 m/s is unaltered. The high energy LLq^i^ 
transition with vf = 1.11 x 10 6 m/s however shows a 
clear change: the structure is much sharper in the second 
measurement. In Fig. [SJa the contribution of the LL_i^o 
transition to <j_(uj) is plotted for both measurements. 
Similar to LLo-vi^, it shows a sharpening in the second 
measurement. The fact that both LLo-j-i^ and LL_i_>.o 
transition peaks are changed with environmental doping 
suggests that these transitions originate from the bands 
in the outmost graphene layer, which is most affected 
by surface contamination. The LLo_>.i ia transition with 
the lowest vp is unchanged by the environmental doping, 
therefore it probably comes from deeper graphene layers. 

Fig. [3^ shows that the spectral weights of the LLo^i^ 
and LL_!_j.o transitions are approximately equal, which 
indicates a balance between electrons and holes in the 
top layer. The Fermi velocity of the electrons (1.11 x 10 6 
m/s) in that layer is slightly larger than the one of the 
holes (1.09 x 10 6 m/s). Similar results were found for 
monolayer exfoliated^! and few layer CVEESI graphene. 
The electron-hole asymmetry in our work (about 2%) 
is much smaller than the one found in Ref. 1361 which 
might be due to different rotation angles between the lay- 
ers in the used samples. Indeed, the asymmetry between 
electrons and holes was experimentally shown to depend 
on the relative rotation between subsequent graphene 
layers^, where a large (small) rotation between the lay- 
ers gives a small (large) asymmetry. 



IV. DISCUSSION 

Multilayer graphene on the C-side of SiC is often re- 
garded as a stack of twisted monolayers, electronically de- 




FIG. 6. Total spectral weight of the LLo->i, a , LLo->i,6 and 
LL-i-s-o transitions (circles) as compared to the theoretical 
prediction for the spectral weight of the sum of the LLo-n 
and LL_n.o transitions in ideal monolayer graphene (solid 
line) with the chemical potential between E_i and Ei. 



coupled from each other due to random rotational stack- 
ing, where the doping level varies across the layers due 
to the effects of the substrate and the surface contamina- 
tion. Our observation of the CR and the LL transitions 
(LLi_>2, LLo->.i, LL-i-^o) in the same spectra is a clear 
indication of this doping variation. Indeed, since only 
transitions between occupied and empty states can oc- 
cur, their simultaneous activation implies different posi- 
tions of the chemical potential with respect to the Dirac 
point energy. However, the activation of different LL 
transitions may also be caused by a spatial doping inho- 
mogencity. 

The nearly perfect square-root dependence of the LL 
tran sitio n energies on magnetic field found in previous 
worlPEH anc f confirmed by the present measurements 
(Fig. [2^,b) is a signature of massless Dirac dispersion 
inherent to monolayer graphene. The s ame field de- 
pendence of LLs is obtained by STM^HH. Accordingly, 
ARPES measurements^ show multiple Dirac cones from 
individual layers. However, a number of our observations 
is difficult to fit into a simple picture of completely iso- 
lated monolayers. One of them is a significant (about 
10%) spread of the Fermi velocity in the same sample. 
The largest values vf = 1.09 - 1.11 x 10 6 m/s we at- 
tribute to the outmost graphene layer, based on the ef- 
fect of the surface contamination, while the smallest val- 
ues, Vp = 1.01 - 1.02 x 10 6 m/s, correspond to the 
inner layers, where the transitions LLo-^a and LLi_>.2 
are active due to a weak electron doping. It was the- 
oretically predicte d 16 ^ 18 * 19 ! that in rotationally stacked 
graphene layers the effect of the interlayer interaction 
is to reduce the Fermi velocity with respect to its 'bare' 
value in monolayer graphene. The variation of vf that we 
find can thus be attributed to the effect of layer twisting 
and a randomness of the rotation angles between various 
layers. 

It is important to compare not only the energies but 
also the optical spectral weig hts W of the LL transitions 
with theoretical expectation d 24 * 34 l Experimentally, we 
obtain the weights from the multi-component spectra fit- 
ting, described in Section III A Let us consider the spec- 
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tral weights of the LLo_>i and LL_i_>o transitions. As- 
suming that the chemical potential is between the first 
electron and hole LLs (E_x and Ei), the total weight 
of these transitions in monolayer graphene is given, ac- 
cording to the Kubo formalism for non-interacting Dirac 
fermions, by the transition energy itselP^: 

hW/a = 2{Ex - E ) = 2^2ehv 2 F \B\ (13) 

This dependence on magnetic field for vp — 1.0 x 10 6 
m/s is plotted in Fig. [6] as a solid line. The symbols 
show the sum of the experimental spectral weights of the 
LLo-i.i, a , LLo->-i.fc and LL_i_>.o transitions. The exper- 
imental points are clearly below the theoretical curve. 
From here it follows that even if only in one graphene 
layer the chemical potential is between E_i and Ei 
then the total spectral weight observed is about 2 times 
smaller than the theoretical expectation. In reality it is 
likely that the number of layers in the present sample 
that satisfy the condition for the chemical potential is 
larger, which would make the deviation even stronger. 

A possibly related experimental observation is the 
presence of the optical absorption background found 
from the fitting results of the phenomenological cyclotron 
multi-component model. This component, with vanish- 
ing w c , has a substantial spectral weight spread over a 
broad frequency range. The existence of the background 
shows that a significant amount of the charge carriers 
in the graphene layers neither fall into well defined CR 
nor LL transitions and signals a departure from the iso- 
lated monolayer description. In terms of the optical sum 
rule, the missing spectral weight of the LL transitions is 
transferred to the background. An intriguing question 
is whether this transfer is caused by interlayer coupling 
or by many-body effects within individual layers, such as 
electron-electron and electron-phonon interactions. 

One should notice that the mobility and density of car- 
riers in epitaxial graphene show a large variation from 
sample to sample, even when they are prepared under 
similar conditions®!. This shows that the interlayer twist 
angle, which is at the moment difficult to control exper- 
imentally, is a crucial parameter affecting the electronic 
and therefore optical properties of epitaxial graphene. 

V. CONCLUSIONS 

Using magneto-optical infrared Hall spectroscopy, the 
charge dynamics in multilayer epitaxial graphene grown 
on the C-side of SiC was studied in magnetic fields up 
to 7 T. The diagonal and the Hall conductivities were 



extracted from the absorption and the Faraday rotation 
spectra, respectively. The latter is sensitive to the sign 
of charge carriers that allowed us to distinguish electrons 
and hole like transitions. The mobility and charge den- 
sity of electrons were found, which makes this technique 
a useful contactless characterization tool. 

A general multi-component model, Eq. fl7b 3 provided 
excellent fits at each field to a xx (oj) and a xy \oj) simulta- 
neously. This analysis revealed the coexistence of optical 
transitions between individual LLs with a square-root de- 
pendence of the transition energies on magnetic field as 
expected for isolated monolayer graphene and a quasi- 
classical CR showing a linear magnetic field dependence. 
This is a clear indication of the doping variation across 
the layers. 

We find the simultaneous presence of at least two dis- 
tinct peaks due to transitions between Landau levels 
and 1 (electrons). The separation between these peaks 
corresponds to a difference between the Fermi velocities 
of about 10%. One more peak due to a transition be- 
tween LLs -1 (holes) and is observed. The effect of the 
surface contamination on the spectra tells that both elec- 
trons and holes are present in the top layer and that elec- 
trons have slightly higher Fermi velocity than the holes 
(by 2%). The variation of the Fermi velocity is proba- 
bly related to random twisting angles between graphene 
layers. 

The spectral weight of the LL transitions is shown to 
be significantly reduced with respect to the theoretical 
expectation for a stack of fully decoupled graphene mono- 
layers, assuming the picture of non-interacting electrons 
within each layer. This missing spectral weight correlates 
with the presence of an unexpected broadband optical 
absorption, which is also inconsistent with this simpli- 
fied theoretical model. 

Although the transition energies of the LL transitions 
clearly follow the magnetic field dependence expected for 
isolated graphene, in order to come to a complete pic- 
ture of the complex electronic structure of multilayer 
graphene we need to understand the variation of the 
Fermi velocity, the small optical spectral weight of the 
LL transition and the broad absorption background in 
relation to the twist angle of the layers. Therefore, a 
systematic study of samples where this angle is experi- 
mentally controlled, is required. 
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ence Foundation (SNSF) by the grants 200021-120347 
and IZ73Z0-128026 (SCOPES program), through the 
National Centre of Competence in Research Materials 
with Novel Electronic Properties-MaNEP. We thank S.G. 
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